AO-AO<M»  850 


UNCLASSIFIED 


aerospace  CORP  el  SE6UN00  CALIF  IVAN  A 6ETTIN6  LABS  F/6  16/2 

A rational  method  for  design  of  feedback  control  to  limit  MISSI— ETC(U) 
AU6  77’  D H PLATUS#  J WOLKOVITCH  F04701-76-C-0077 


TR-0077(2240)-4 


SAMSO-TR-77-193 


REPORT  SAMSO-TR-77-J93 


o 

lO 

00 

o 

< 

Q 


Cu 

o 

CJ> 


C-D 

ra  Cl::;::* 
*a!s:  c:::^ 


A Rational  Method  for  Design  of  Feedback 
Control  to  Limit  Missile  Dispersion 


Aerophysics  Laboratory 
The  Ivan  A.  Getting  Laboratories 
The  Aerospace  Corporation 
El  Segundo,  Calif.  90245 


15  August  1977 


Interim  Report 


Prepared  for 


SPACE  AND  MISSILE  SYSTEMS  ORGANIZATION 
AIR  FORCE  SYSTEMS  COMMAND 
Los  Angeles  Air  Force  Station 
P.O.  Box  92960,  Worldway  Postal  Center 
Los  Angeles,  Calif.  90009 


r 


rn  rp 

OCT  4 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED 


This  interim  report  was  submitted  by  The  Aerospace  Corporation, 

El  Segundo,  CA  90245,  under  Contract  No.  F0470 1-76-C-0077  with  the  Space 
and  Missile  Systems  Organization,  Deputy  for  Advanced  Space  Programs, 
P.O.  Box  92960,  Worldway  Postal  Center,  Los  Angeles,  CA  90009.  It  was 
reviewed  and  approved  for  The  Aerospace  Corporation  by  W.  R.  Warren,  Jr. 
Director,  Aerophysics  Laboratory.  Lieutenant  Batki,  SAMSO/YAPT,  was 
the  project  officer  for  Advanced  Space  Programs. 

This  report  has  been  reviewed  by  the  Information  Office  (OI)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS, 
it  will  be  available  to  the  general  public,  including  foreign  nations. 

The  technical  report  has  been  reviewed  and  is  approved  for  publication. 
Publication  of  this  report  does  not  constitute  Air  Force  Approval  of  the 
report's  findings  or  conclusions.  It  is  published  only  for  the  exchange  and 
stimulation  of  ideas. 


Dara  Batki,  Lieutenant,  USAF 
Project  Officer 


Robert  W.  Lindemuth,  Lt  Col,  USAF 
Chief,  Technology  Plans  Division 


FOR  THE  COMMANDER 


Deputy  for  Advanced  Space  Programs 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wh»n  Dmim  Entartd) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.  GOVT  ACCESSION  NO.  3.  RECIPIENT’S  CATALOG  NUMBER 


A RATIONAL  METHOD  FOR  DESIGN  OF 
FEEDBACK  CONTROL  TO  LIMIT  MISSILE 
DISPERSION  , 


7.  AUTHORfa; 

Daniel  H./piatus 


Julian/Wolkovitch 


5.  TYPE  OF  REPORT  ft  PERIOD  COVERED 

(,  Inte  rim  j 


Tr''-''(^ 7 7 ( lZ4^-°4  ; 


8.  CONTRACT  OR  GRANT  fJUHBERCa; 

/’"1fM70  1 - 76- C- 0077 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

i/^he  Aerospace  Corporation 
El  Segundo,  Calif.  90245 


10.  PROGRAM  element.  PROJECT,  TASK 
AREA  8 WORK  UNIT  NUMBERS 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


12.  REPORT  DATE 


Space  and  Missile  Systems  Organization  / j ) 

Air  Force  Systems  Command 

Los  Angeles,  Calif.  90009  


14  monitoring  agency  NAME  4 ADDRESSflf  di//or«n(  Irom  Controlling  Otlico) 


16.  DISTRIBUTION  ST  ATEMEN  T fof  (fila  ReporlJ 


15  Aug^  mn 


1 5a.  DECLASSIFICATION  ^DOWNGRADING 

schedule 


Approved  for  public  release;  distributi  .mited 


17.  DISTRIBUTION  STATEMENT  (ol  the  obBlracI  aniarad  In  Block  20,  II  dlllaranl  horn  Report) 


KEY  WORDS  (Continue  on  reverse  side  if  necessery  and  Identify  by  block  number) 

Missile  Control 
.Missile  Accuracy 
Missile  Dispersion 
Missile  Aerodynamics 


.^ABSTRACT  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 

" A method  for  the  design  of  feedback  control  to  limit  cross-range  dispersion  of 
spinning  missiles  is  described.  The  missile  transverse  velocity  in  the  cross 
plane  is  selected  as  the  cost  function,  and  optimal  linear  feedbacks  are  de- 
rived that  minimize  this  velocity  in  the  presence  of  aerodynamic  disturbances 
that  cause  lift  nonaveraging.  A closed-form  solution  is  obtained  for  the  open- 
and  closed-loop  missile  response  to  simple  disturbance  moments.  Cross- 
range dispersion  control  complements  control  of  range  error  by  >7  .9 


• FACSIMILE! 


UNCLASSIFIED  /yLO*' 

security  CLASSIFICATION  OF  THIS  PAGE  (IWi»n  i 


UNCLASSIFIED 


security  classification  of  this  PAOerlWl«n  Dmim  Bnt»nd) 


security  classification  of  this  RAOEflThMi  DMA  Enimnd) 


PREFACE 


The  authors  are  grateful  to  M.  E.  Brennan 
numerical  computations. 


for  her  assistance  with  the 


I 


CONTENTS 


PREFACE  1 

I.  INTRODUCTION 7 

II.  ANALYSIS 9 

A.  Control  Equations 9 

B.  Solution  for  Linear  Optimal  Control  14 

C.  Response  to  Impulsive  Yaw  Moment 17 

III.  NUMERICAL  EVALUATION 21 

IV.  SYSTEM  IMPLEMENTATION 29 

V.  CONCLUSIONS 31 

REFERENCES  33 

NOMENCLATURE 35 


-3- 

L-  I 


FIGURES 


1.  Transverse  velocity 10 

2.  Open-loop  angle-of-attack  response  to  impulse 22 

3.  Open-loop  precession  rate  response  to  impulse 23 

4.  Open-loop  dispersion  velocity  response  to  impulse 24 

5.  Closed-loop  angle-of-attack  response  to  impulse  25 

6.  Closed-loop  precession  rate  response  to  impulse 26 

7.  Closed-loop  dispersion  velocity  response  to  impulse 27 

TABLE 

1.  Net  dispersion  velocity  in  ft/sec 28 


-5- 


I.  INTRODUCTION 


Dispersion  of  missiles  can  be  separated  into  two  principal  categories; 

(1)  cross-range  dispersion  resulting  from  lift  nonaveraging  and  (2)  up-  or 

down-range  error  associated  with  drag  uncertainty  and  atmospheric  variations 

such  as  winds.  Winds  also  contribute  to  cross-range  dispersion.  Control 

of  up-  or  down-range  error  by  compensating  for  drag  uncertainty  (drag 

modulation)  has  been  treated  previously.  It  has  been  shown  that  very  little 

energy  is  required  to  effect  large  changes  in  drag  by  yaw-moment  control 
1 2 

of  angle  of  attack.  ’ The  missile  is  controlled  in  a circular  coning  motion 
in  which  the  drag  is  a strong  function  of  the  coning  half-angle  (angle  of 
attack).  Changes  in  drag  caused  by  extraneous  sources,  as  detected  by  drag 
deceleration,  can  be  compensated  for  by  appropriately  increasing  or 
decreasing  the  angle  of  attack. 

Cross- range  dispersion  resulting  from  lift  nonavoraging  is  one  of  the 

greatest  contributors  to  missile  impact  error.  The  sources  of  such  dispersion 

are  quite  varied  and  include  launch  errors  caused  by  muzzle  disturbance,  ^ 

4 5 6 7 

small  mass  and  configurational  asymmetries,  ’ ’ ’ and  asymmetric 

g 

boundary-layer  transition  (in  the  case  of  ballistic  reentry  vehicles).  It  is 

well  known  that  dispersion  caused  from  variation  in  body-fixed  asym- 

5 8 

metrics  varies  approximately  inversely  with  roll  rate,  ’ and  that  such  dis- 
persion can  be  controlled  by  maintaining  the  roll  rate  at  a sufficiently  large 
steady  value  (roll  control).  However,  certain  forms  of  distvirbance  moments 
can  adversely  affect  the  lift-vector  precession  rate  and  produce  lift  non- 

g 

averaging  dispersion  somewhat  independent  of  the  roll  rate. 

In  this  report,  a general  method  is  described  for  the  design  of  feedback 
control  to  limit  cross- range  dispersion  caused  by  lift  nonaveraging.  This 
control  would  complement  drag  control  to  limit  up-  and  down- range  error: 
the  resulting  two-loop  system  would  control  major  sources  of  missile 
dispersion. 
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II.  ANALYSIS 


A.  CONTROL  EQUATIONS 

Cross- range  dispersion  resulting  from  lift  nonaveraging  can  be 
described  in  terms  of  the  complex  missile  transverse  velocity  V = v + iw  in  a 
plane  normal  to  the  average  flight  path.  This  velocity  is  defined  in  terms 
of  the  total  angle  of  attack  0 and  lift- vector  precession  angle  4j  according 
to^  (Fig.  1). 

v(t)  = v(0)  - ee^'^dt  (1) 

0 

where  Lgis  the  lift  force  derivative  and  m is  the  missile  mass.  If  the 
avc  '■  value  of  the  transverse  velocity  V(t)  is  set  equal  to  zero  prior  to  any 
■nee  in  9 and  4j  that  causes  lift  nonaveraging,  then  the  cross- range 
,n  resulting  from  a disturbance  is  proportional  to  the  magnitude  of 
.average  value  of  V at  sufficiently  large  time  t after  the  disturbance. 

The  behavior  of  the  angles  6 and  is  described  approximately  by  the 
undamped  equations  of  motion  for  a slowly  rolling  (p  — const.),  axisymmetric 
missile  ^ 


e + (uo‘ 


• 2 

4^  ) e = ^ 


M 


6p 


(2) 


M M 

+ 294j  = + —4^ 


(3) 
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where  M and  M are  aerodynamic  pitch  and  yaw  disturbance  moments  in 

P y * 

the  wind- referenced  axes  and  and  are  applied  control  moments. 

The  control  problem  is  to  define  the  moments  Mg^  and  Mg^  in  terms  of 
appropriate  state  variables  that  minimize  the  net  transverse  velocity  V 
resulting  from  some  disturbance.  The  problem  is  that  of  an  optimal  regulator 
with  a cost  function  equal  to  the  magnitude  of  the  integral  in  Eq.  (1). 

It  is  assumed  that  the  missile  is  initially  untrimmed  in  a circular 
coning  motion  to  permit  a sufficient  level  of  drag  control  by  either  increasing 
or  decreasing  the  angle  of  attack.  ^ Cross-range  dispersion  is  then  con- 
trolled by  minimizing  the  coupling  between  angle-of-attack  and  precession 
angle  perturbations  that  cause  lift  nonaveraging  in  accordance  with  Eq.  (1). 
The  equations  of  motion,  Eqs.  (2)  and  (3),  can  be  linearized  in  terrns  of 
small  perturbations  about  the  quasi- steady  values  9=9  and  = u)^.  With 
0 =9+9^  and  4^  = + ‘^^4.  substituted  in  Eqs.  (2)  and  (3)  and  higher-order 

terms  neglected  (where  0+  and  4^+  are  small  perturbations),  the  resulting 
linear  control  equations  can  be  written 

M 

- 2w  04j  + = -^  + (4) 

+ n ^ I I 


2co  e ^ +0^^ 
n + + 


M 

y 

I 


M 


6y 


(5) 


where  the  control  moments  are  assumed  to  be  linear  functions  of  the  state 
variables  9^,  6^,  and  4^^,  according  to 


Aerodynamic  damping  is  omitted,  and  the  roll  rate  is  assumed  to  be  small 

such  that  I p/2I  « w . 

n 


-11- 


(6) 


M 


= - AG^  - BG^  - 


M 


5y 


DG EG , 

T + + 


(7) 


It  is  convenient  to  nondimensionalize  the  control  equatic  ».s  by  defining  new 
variable  s 


Q'=0,/G,  \ =4',/w  , T =00  t 
+ + n n 


(8) 


If  these  variables  are  subsitituted  in  Eqs.  (4)  through  (7),  sinc'-  d/dt  = oo^d/dr, 
and  the  Laplace  transform  with  respect  to  t defined  by 


- ST 

( ) e dr 


.n-l 

is  taken,  the  control  equations  can  be  written 


■ 2 

s +bs+a  c-2 

Of 

m 

r: 

P 

(2  + e)s+d  s+f 

\ 

m 

L yj 

where  the  nondimensional  feedbacks  a,  b,  c,  ...  are  defined  by 


A.B  C.D  E.F 

a s — =.  b = — , c = =- , d s c s — , f = — ^ 

^ 00  oofi  2 OO  OO0 

oo  n n 00  n n^ 

n n 


(9) 


(10) 


(11) 


This  step  is  not  essential  for  the  analysis,  but  it  facilitates  the  numerical 
evaluation  described  later. 
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and  the  nondimensional  transformed  disturbance  moments  are 


m = 
P 


w^ei 

n 


^ 7'{M  ),  m = — 

p V " 


w GI 
n 


(12) 


The  transformed  nondimensional  control  moments,  Eqs.  (6)  and  (7),  can 
similarly  be  written 


mj.  = -[  (a  + bs)  a (s)  + c\(s)]  (13) 

6p 


m,  = -[  (d  + es)  a (s)  + fX(s)]  (14) 

6y 

The  transverse  velocity,  Eq.  (1),  can  be  written  in  nondimensional 
form  with  the  substitutions  G = G(1  + »)  and 


(oj  + ^ ,)  dt 
n + 


1 + X)  dr 


(15) 


which  yields 


V(t)  = V(0)  - if 


w er' 

n 


( 1 + O' ) exp 


r*/ 

L-'o 


(1  + X)dT 


dr 


(16) 


2 . 

where  the  relation  for  pitch  frequency  = L^x^^/I  is  used  and  a characteristic 
length  I = I/mx^^  is  defined.  The  control  problem  defined  by  Eqs,  (10)  and 
(16)  is  that  of  a linear  optimal  regulator.  The  optimal  linear  feedbacks  a, 
b,  c,  ...  that  will  minimize  the  cost  function,  Eq.  (16),  are  to  be  deter- 
mined, It  is  necessary  to  include  the  additional  constraint  on  the  cost  func- 
tion that  the  control  moments  remain  within  acceptable  limits  to  preclude  the 
trivial  case  of  very  large  feedbacks  by  redefining  the  cost  function 
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I = WjVV="  + W2 


I 2 2, 


where  is  the  complex  conjugate  of  V,  and  are  the  nondimensional 

control  moments  defined  by  Eqs.  (13)  and  (14),  and  , W^.  and  are 
suitable  weighting  functions. 

B.  SOLUTION  FOR  LINEAR  OPTIMAL  CONTROL 

The  solution  to  Eq.  (10)  for  a(s)  and  X(s)  is  of  the  form 


N,(s) 


N2(s) 


“"nriF  ’ D(s) 


where  N^(s),  N2(s),  and  D(s)  are  polynomials  in  s of  order  depending  on  the 
form  of  the  disturbance  moments  m^  and  m^  and  with  coefficients  containing 
the  feedbacks  a,  b,  c,  ...  to  be  determined.  The  control  moment  terms  in 
Eq.  (17),  in  the  form  of  integral- square  values,  can  be  readily  evaluated 
from  the  results  of  Eq.  (18)  in  terms  of  Phillips  integrals. 

For  example,  if  the  system  is  third  order  and  the  pitch-control 

moment  m,  , Eq.  (13),  has  the  form 
6p 


c^s  + c ^s  + Cq 


d s^  + d^s^  + djS  + dj. 


/°°  2 

m^pdr  has  the  value 


c^  d^dj  + (c^  - 2cQC2)dqd3  + Cq  d2d3 

2d.d„(d.d'  - d„d„) 


^0^3'  1 2 
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Evaluation  o£  the  transverse  velocity  term  VV="  is  much  more  difficult 
because  of  the  nonlinear  coupling  between  a and  \ in  Eq.  (16).  A reasonable 
approximation  to  Eq.  (16)  can  be  obtained,  however,  by  writing  it  in  the  form 


V(t)  = V(0)  - if  cj^e [ 1 + o(t)]  e’'^l^^^e^^dT 


where 


\j(r)  =y  \(T)dr,  \^{s) 


and  the  upper  limit  of  the  integral  in  Eq.  (21)  is  taken  to  be  sufficiently  large 
to  include  the  net  change  in  V(t)  resulting  from  perturbations  a (t)  and  K(t). 
Without  loss  of  generality,  this  upper  limit  can  be  taken  as  co,  and  Eq.  (21) 
becomes  the  Fourier  transform  of  the  function  f(r) 


V = V(0)  - if  w qT  f(r)  e^^dr 


whe  re 


f(T)  ::  [1  +a(T)]e'^l^^^  (24) 

Since  o'(t)  and  \(t)  are,  in  general,  small  pe rtu rbations , according  to  their 
definitions,  Eq.  (8),  the  exponential  in  Eq.  (24)  can  be  expanded  to  obtain  for 
f(r)  the  app  rerXimation 

1 2 

f(T)  « 1 +a(r)  + i\j(r)  + io(T)\j(T)  _ _ [ \j(t)]  + . . . (25) 


Equation  (23)  can  be  evaulated  in  terms  of  the  Laplace  transform  of  f(T), 

Eq.  (9),  with  s = -i.  The  constant  term  in  Eq.  (25)  is  the  contribution  to 

V(t)  from  6 prior  to  a disturbance  that  causes  ©(t)  and  X(t).  If  V(0)  in  Eq.  (23) 

is  defined  as  the  average  value  of  V(t)  prior  to  a disturbance,  both  the 

constant  and  V(0)  can  be  dropped,  since  only  the  net  change  in  V from  the 

8 

perturbation  is  of  interest.  Equation  (23)  can  then  be  written 


AV  = -i/cj  6'/'  o(t)  + i\  (t)  + io'(T)  X. . (t)  - 4-[  \ . (t) ] ^ + 
n 1 1 1 


s = -1 


-i?  to  9 
n 


a{s)  +x^  r i//|o(T)  \^(T)]  - J (T)]^  + 


S = -1 


(2  6) 


where  O'(s),  \{s) , and  X-j(s)  are  given  by  Eqs.  (18)  and  (22).  A first-order 
solution  consists  simply  of  the  first  two  terms,  which,  from  Eq.  (18),  can 
be  written 


N 


AV  = -ifw  e 

1st  ord 


(-i)  - N (-i) 

dTI) 


(27) 


A second-order  solution  includes,  in  addition  to  Eq.  (27),  the  terms  involving 
the  Laplace  transform  of  a product.  For  the  case  where  o(s)  and  X(s),  given 
by  Eqs.  (18)  and  (22),  consist  of  m and  n first-order  poles,  respectively, 
and  no  others,  the  Laplace  transforms  are  given  by^^ 


'/'[  o(T)  Xj(T)] 


m 


(s  - s^) 


(28) 


2 ^ ^2^"=^ 


j = l V j 


(29) 
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where  s,  are  roots  of  D(s)  = 0 and  s.  are  roots  of  D.(s)=  sD(s)  - 0.  Thus, 
k J 1 

the  approximate  solution  for  V,  given  by  Eqs.  (26)  through  (29),  is  obtained 
directly  from  the  s-plane  solution  of  the  control  equations  in  the  form  of 
Eqs,  (18)  and  (22).  The  optimum  feedbacks  are  those  values  that  minimize 
the  cost  function,  Eq.  (17),  according  to 


II 

9 a 


= 0, 


(30) 


C.  RESPONSE  TO  IMPULSIVE  YAW  MOMENT 

The  method  is  illustrated  with  a simple  example  that  is  amenable  to  an 
analytical  solution.  The  response  to  an  impulsive  yaw  monient 


M (t)  = M6(t)  (31) 

y 


is  considered  such  that  m = M/w  Si,  from  Eq.  (12),  is  a constant  and 

t y ^ 

M = 0.'  The  solution  to  Eq.  (10)  for  »(s)  and  \(s)  in  the  form  of  Eq.  (18) 
P 
is 


o(s) 


m (2  - c) 

y 

D(s) 


(32) 


\(s) 


rn  ( s 

y 


2 


+ bs 

WT 


+ a) 


(33) 


where 


D(s)  = s^  + (b  + f)s^  + [ a + bf  + (2  + e)(2  - c)]  s + af  + d(2  - c) 


(34) 


^Since  the  unit  impulse  function  5(t)  has  units  of  1/t,  it  must  be  nondimension- 
alized  by  dividing  by  oj  prior  to  h>king  the  _Laplace  transform  defined  in  Eq. 
Eq. 


(9),  which  accounts  for  the  form  M/w^0I 


rather  than  M/w^  61  for  m 
n y 


J 
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The  Routh  stability  criterion^^  requires  that  the  feedbacks  satisfy  the 
relations 


af+d(Z-c)>0 
a + bf  + (2  + e)  (2  - c)  >0 
b + f >0 

(b  + f)  [ bf  + (2  + e)  (2  - c)l  + ab  - d(2  - c)  >0  (35) 


1 1 


! 

These  are  satisfied  with  the  feedbacks  a,  f ^ 0 -..nd  b = c = d=  e=0.  Optimum 
I values  of  these  two  feedbacks  will  be  found  that  minimize  dispersion  in 

i accordance  with  our  cost  function,  Eq.  (17).  The  perturbation  solution, 

j Eqs.  (32)  through  (34),  reduces  to 

[ 

1 

' 2m  m (s^  + a)  _ _ 

a(s)^-y,  \(s)  = — X . D(s)  = s + fs  + (a  + 4)  s + af  (36) 

D(s) 

and  the  first-order  solution,  Eq.  (27),  is  found  to  be 


AV 


1st  o rd 


(Mf  /I)(3  - a)  f 3 + a - i(a  - l)f] 

2 2 2 
(a  - 1)  f + (a  + 3) 


(37) 


This  expression  vanishes  for  a = 3,  independent  of  f.  Hence,  to  the  first- 
order  approximation,  dispersion  caused  from  an  impulsive  yaw  moment  can 
bo  eliminated  with  two  simple  feedbacks.  In  a more  general  case,  the 
feedbacks  would  be  determined  to  minimize  both  the  dispersion  term  VV='= 
and  the  control  moments  in  our  cost  function,  Eq.  (17),  with  suitable 
weighting  of  each  term.  Since  the  feedback  a = 3 eliminates  dispersion,  to 
the  first-order  approximation,  and  a nonzero  value  of  the  feedback  f is 
required  for  stability,  the  magnitude  of  f is  determined  to  effect  equal 


j 
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integ ral- square  values  of  the  control  moments  in  order  to  illustrate  the 
method.  The  control  moments,  Eqs.  (13)  and  (14),  with  a = 3,  b = c=  d = 
e = 0,  and  the  solution,  Eq.  (36),  can  be  written 


2a  m 
^6p  D(s) 


y 


6y 


(s^  + a) 

D(i) 


(38) 


where  D(s)  is  defined  in  Eq.  (36).  The  integral- square  values  of  these 

moments  are  given  by  the  integral,  Eq.  (20),  where  the  coefficients  c^  and 

d are  defined  for  each  moment  by  Eqs.  (36)  and  (38).  If  the  numerators  of 
n 

Eq.  (20)  are  equated  for  the  two  moments  (since  the  denominators  are  the 
same),  the  nondimensional  feedbacks  a and  f must  satisfy  the  relation  f = \Ta 
for  equal  integ  ral- square  values  of  the  control  moments.  Thus,  the  optimal 
linear  first-order  values  of  our  two  feedbacks  for  the  assumed  cost  function 
are  a = 3 and  f = \T3.  These  values  should  result  in  zero  dispersion,  to  the 
first-order  approximation,  and  require  equal  integral- square  values  of  the 
control  moments.  The  open-loop  dispersion  caused  by  an  impulsive  yaw 
moment  M is  readily  obtained  from  Eq.  (26)  and  the  perturbation  solution, 
Eqs.  (32)  through  (34),  with  all  of  the  feedbacks  equal  to  zero.  The  first- 
order  approximation  to  this  dispersion  is  found  from  Eq.  (37)  with  a = f = 0 
to  give  the  simple  result 
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III.  NUMERICAL  EVALUATION 

The  equations  of  motion,  Eqs.  (2)  and  (3),  were  integrated  numerically 
to  obtain  both  the  open-  and  closed-loop  responses  to  an  impulsive  negative 
yaw  moment  with  the  simple  optimal  feedbacks  derived  above.  The  results 
are  compared  with  first-  and  second-order  analytical  approximations  to  the 
dispersion  velocity.  Angle  of  attack,  precession  rate,  and  transverse 
velocity  histories  are  shown  in  Figs.  2 through  7 for  the  open-  and  closed- 
loop  responses.  The  inputs  and  system  parameters  are  as  follows: 

I =20  slug  ft^ 

f.  = 6.  566  ft 

CO  =40  rad /sec 

_n 

0 =2  deg 

a =3 

f = ^f3 

b =c  = d = e = 0 

M = /*  M dt  = -6.092  ft  lb  sec 

The  feedbacks  used  in  the  simulation  of  Figs.  5 through  7 (a  = 3,  f = \/  3)  were 
determined  to  give  zero  net  dispersion  velocity  on  the  basis  of  the  first-order 
approximation,  Eq.  (37).  The  closed-loop  velocity,  computed  from  a 
numerical  integration  of  the  nonlinear  equations  (2)  and  (3),  is  shown  in 
Fig.  7.  The  feedback  control  limits  the  dispersion  to  approximately  3%  of 
the  open-loop  value.  The  second-order  contributions  to  the  net  dispersion 
velocity,  Eq.  (26),  can  be  computed  from  the  relations,  Eqs.  (28)  and  (29), 
for  the  Laplace  transform  of  a product,  and  from  the  control  solutions, 

Eq.  (36).  The  second-order  terms  are  found  to  be 
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VELOCITY  (w),  ft/sec 


VELOCITY  (w), ft/sec 


-8-4048 
VELOCITY  (v),  ft/sec 

Fig.  7.  Closed-loop  dispersion  velocity  response  to  impulse 
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= 0.0222  - 0.01231  ft/sec 


(41) 


which  should  be  added  to  the  first-order  solutions,  Eqs.  (39)  and  (37),  of 
-2.00  ft/scc  and  zero,  respectively.  The  results  are  summarized  in 
Table  1,  which  is  a comparison  of  the  open-  and  closed-loop  theoretical 
values  with  those  computed  from  a numerical  integration  of  the  nonlinear 
equations  of  motion. 


Table  1.  Net  dispersion  velocity  in  ft/sec 


Open 

Closed 

Loop 

Loop 

First-Order  Approximation 

2.00 

0 

First-  and  Second-Order 
Approximations 

1.97 

0.025 

Numerical  Integration 

2.00 

0.067 
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IV.  SYSTEM  IMPLEMENTATION 


Implementation  of  a control  system  to  sense  the  required  state 

variables  and  to  generate  wind- referenced  pitch  and  yaw  control  moments 

1,2 

has  been  discussed  in  conjunction  with  drag  control  systems.  ’ The  princi- 
pal feature  of  such  systems,  which  differs  from  that  of  more  conventional 
aerodynamic  control  systems,  is  the  requirement  to  generate  aerodynamic 
control  moments  in  the  wind  plane,  which  rotates  relative  to  body-fixed 
axes  at  a rate  approximately  equal  to  the  difference  between  the  roll  rate  and 
the  natural  pitch  frequency  when  the  vehicle  is  untrimmed.  Because  relatively 
small  control  moments  are  required  to  correct  for  aerodynamic  disturbances 
that  would  otherwise  cause  lift  nonaveraging,  the  power  requirements  are 
small  enough  that  wind-referenced  moments  can  be  generated  by  modulation 
of  body-fixed  control  surfaces  or  reaction  jets  at  the  windward-meridian 
rotation  frequency.  The  state  variables  selected  for  this  analysis  consist 
of  the  angle  of  attack,  the  pitch  rate,  and  the  precession  rate.  Since  the 
transverse  dispersion  velocity,  Eq.  (1),  is  the  integral  of  the  lateral 
acceleration,  a more  appropriate  control  variable  than  total  angle  of  attack 
is  lateral  acceleration,  which  is  measured  directly  from  a resolution  of 
body-fixed  lateral  accelerometers.  The  pitch  rate  6 can  similarly  be 
obtained  from  a resolution  of  body-fixed  lateral  rate  gyro  measurements. 

The  precession  rate  is  the  difference  between  the  roll  rate  and  the  windward- 
meridian  rotation  frequency.  With  the  missile  untrimmed  in  a circular  coning 
motion  at  a quasi-steady  angle  of  attack  6 , a body-fixed  lateral  accelero- 
meter will  measure  a strong  periodic  signal  at  the  windward-meridian 
rotation  frequency.  The  difference  between  this  signal  frequency  and  the 
roll  rate  is  the  precession  rate.  Hence,  all  required  control  parameters 
can  be  obtained  from  conventional  strapped-down  sensors. 
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For  control  of  yaw-moment  disturbances  only,  as  discussed  in  the  numerical 
example,  the  only  feedbacks  required  are  angle  of  attack  (lateral  accelera- 
tion) and  precession  rate.  For  a small  or  known  roll  rate,  control  could  be 
achieved  with  information  derived  solely  from  lateral  accelerometel  2. 


V.  CONCLUSIONS 


Cross- range  dispersion  of  spinning  missiles  resulting  from  lift 
nonaveraging  can  be  controlled  by  application  of  wind-referenced  pitch  and 
yaw  moments,  A linear  optimal  control  is  derived  that  minimizes  transverse 
dispersion  velocity  selected  as  the  cost  function.  An  approximation  to  the 
dispersion  velocity  pennits  the  optimal  feedbacks  to  be  obtained  directly 
from  the  s-planc  solution  of  linearized  control  equations  by  means  of 
classical  techniques.  The  method  also  permits  closed-form  solutions  for 
the  open-  and  closed-loop  system  response  to  simple  disturbances.  The 
control  loop,  defined  to  minimize  cross-range  dispersion,  complements 
control  of  up-  and  down-range  error  by  angle-of-attack  control  of  drag.  The 
combined  two-loop  system  can  limit  major  sources  of  missile  dispersion. 


L . 
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po|rfi»nnMa  paGB  MOft  IXXMED 
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NOMENCLATURE 


A,  B,  C,  D,  E,  F 


feedback  gains 


a,  b,  c,  d,  e,  f 


transformed  feedback  gains,  Eq,  (11) 


"^l'  ^2’  *^0’  ’^2 


m , m 
p y 


""6p'  ”"6y 


M , M 

p y 

^6p'  ^6y 


Eq.  (19) 

pitch  or  yaw  moment  of  inertia;  cost 
function,  Eq.  (17) 

roll  moment  of  inertia 

characteristic  length,  I/mx^^. 

lift  force 

lift  force  derivative 
missile  mass 

transformed  pitch  and  yaw  disturbance 
moments,  (Eq.  (12) 

transformed  pitch  and  yaw  control 
moments,  Eqs.  (13)  and  (14) 

pitch  and  yaw  disturbance  moments 

pitch  and  yaw  control  moments 

roll  rate 

Laplace  trans^'orm  variable 
time 

Y-component  of  transverse  velocity 

transverse  velocity,  v + iw 

net  change  in  average  transverse 
velocity 

Z-component  of  transverse  velocity 


Wj,  W3 


weighting  constants 
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st 

Y.  Z 
a 

5(t) 

0 

6 

\ 


4^ 


static  margin  (distance  of  center  of 
pressure  aft  of  center  of  mass) 

cross-plane  coordinates 


e ^/e 


unit  impulse  function  | 

t| 

angle  of  attack  (Euler  angle)  ij 

mean  angle  of  attack 


angle  of  attack  perturbation 


^ + /“n 

Eq.  (22) 


nondimensional  time,  co  t 

n 

precession  angle  (Euler  angle) 
precession  rate 
precession- rate  perturbation 
undamped  natural  pitch  frequency 


I 
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THE  IVAN  A.  GETTING  LABORATORIES 

The  r.aboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems,  Ver* 
satility  and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 

fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory;  Atmospheric  reactions  and  atmos- 
pheric optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory;  Development  of  new  materials:  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory;  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmos^ereV  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 

THE  AEROSPACE  CORPORATION 
El  Segundo,  California 
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